Evidence of Doping-Dependent Pairing Symmetry in Cuprate Superconductors 
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Scanning tunneling spectroscopy (STS) studies reveal long-range spatial homogeneity and pre- 
dominantly d ;c 2_j / 2-pairing spectral characteristics in under- and optimally doped YBa2Cu307_5 
superconductors, whereas STS on YBa2(Cuo.9934Zno.oo26Mgo. 004)306.9 exhibits microscopic spatial 
modulations and strong scattering near the Zn or Mg impurity sites, together with global suppres- 
sion of the pairing potential. In contrast, in overdoped (Yo.7Cao.3)Ba2Cu307_i , (d x 2_ y 2 + s)-pairing 
symmetry is found suggesting significant changes in the superconducting ground-state at a critical 
doping value. 



One of the most important issues associated with the 
cuprate superconductors is the possible existence of a 
quantum critical point (QCP) near the optimal hole dop- 
ing (p « 0.16 per CUO2) in these doped Mott antifcr- 
romagnets |l], |^, [|. Although non-monotonic doping- 
dependent physical properties in the underdoped (p < 
Po) and overdoped (p > p Q ) regimes are suggestive of 
the presence of a QCP, an unambiguous proof would re- 
quire identifying the relevant broken symmetry associ- 
ated with a critical doping level p c . Various theoreti- 
cal studies 0, |J have investigated the dependence of the 
ground state of doped Mott antiferromagnets on the dop- 
ing level and the strength of exchange and Coulomb in- 
teractions 0. The relevant symmetries associated with 
the competing orders include: electromagnetic U(l), spin 
rotation-invariance SU(2), square lattice space group C± v 
and time-reversal symmetry T ^ . Experimentally, while 
the d x 2 _y2 pairing symmetry is known to dominate in the 
optimally doped cuprates []|, ^, 0, |§|, || , possible doping- 
dependent pairing symmetry has been suggested Jl(| . In 
general, the doping dependence of the pairing symmetry 
and the issue of quantum criticality must be considered 
under the premise of spatial homogeneity in the pairing 
potential. Recent finding of nano-scale variations in the 
measured energy gap of Bi 2 Sr 2 CaCu208+ x (Bi-2212) Q 
and the implication of macroscopic clusters in overdoped 
La2- x Sr x Cu04 JT2[ present possible complications in un- 
raveling the existence of a QCP. 

A possible consequence of a QCP is the doping- 
dependent pseudogap phenomenon (ll| , which may rep- 
resent a precursor for superconductivity in the cuprates 
p3[ . Early experiments on the Bi-2212 system reported 
a measured energy gap A*(p) that increased monoton- 
ically with decreasing p and was nearly independent of 
temperature @ [0| |§ [HJ fl|. However, the low- 
temperature spectra of the optimally doped and under- 
doped Bi-2212 appeared to consist of a sharp peak feature 



on top of a broad "hump" . Recent bulk measurements 
on Bi-2212 mesas |^] demonstrated strong temperature 
dependence associated with the sharp peak, which van- 
ished at the superconducting transition temperature T c , 
while the hump feature persisted well above T c . The co- 
existence of these two gap-like features in the supercon- 
ducting state has been attributed to different physical 
origin associated with each gap pT 
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In this work, we address some of these issues via stud- 
ies of the directional and spatially resolved quasiparticlc 
tunneling spectra on the YBa2Cu30y_5 (YBCO) with a 
range of doping levels. The doping dependence of the 
pairing symmetry, pairing potential and spatial homo- 
geneity is derived from these studies. 

The samples used in this investigation included three 
optimally doped YBCO single crystals with T c — 92.9 ± 
0.1 K, three underdoped YBCO single crystals with T c = 
60. Oil. 5 K, one underdoped YBCO c-axis film with T c = 
85.0 ± 1.0 K, two overdoped (Y .7Cao.3)Ba2Cu307_,5 
(Ca-YBCO) c-axis films || with T c = 78.0 ± 2.0 
K, and one optimally doped single crystal contain- 
ing small concentrations of non-magnetic impurities, 
YBa2(Cuo.9934Zno.oo26Mgo.oo4)306.9 ((Zn,Mg)-YBCO), 
with T c = 82.0 ± 1.5 K §, ||. The spectra of YBCO 
single crystals were taken primarily with the quasiparti- 
cles tunneling along three axes: the antinode axes {100} 
or {010}, the nodal axis {110}, and the c-axis {001}; 
while those of the pure and Ca-YBCO films were taken 
along the c-axis. All samples except (Zn,Mg)-YBCO 
are twinned. The surface was prepared by chemical 
etching^, ^4j, and samples were kept either in high- 
purity helium gas or under high vacuum at all times. 
Our surface preparation has the advantage of terminat- 
ing the YBCO top surface at the Cu02 plane by chemi- 
cally passivating the layer while retaining the bulk prop- 
erties of the constituent elements [^3[ thus yielding 
reproducible spectra for samples of the same bulk stoi- 
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chiometry, although direct constant-current mode atomic 
imaging of the chemically inert surface becomes difficult. 
For comparison, surfaces of vacuum-cleaved YBCO sam- 
ples are found to terminate at the CuO-chain layer, which 
are prone to loss of oxygen and development of surface 
states |2^, [24). For vacuum-cleaved Bi-2212, the surface 
typically terminates at the BiO layer. 

Figure illustrates representative tunneling conduc- 
tance (dlNs/dV) vs. voltage (V) data with high spa- 
tial resolution for YBCO samples at 4.2 K: (a) optimally 
doped YBCO, with the average quasiparticlc momentum 
k || {110} and the tip scanning along {001}; (b) un- 
derdoped YBCO, with k || {100} and scanning along 
{001}; and (c) Ca-YBCO, with k || {001} and scanning 
along {100}. The normalized (dlNs/dV) vs. V spectra 
of under- and (Zn,Mg)-doped YBCO for k || {100} are 
shown in the inset of Fig. 0(b), and those of underdoped, 
(Zn,Mg)-doped, and Ca-doped YBCO for k || {001} 
are illustrated in Fig. 0(a). To verify the spatial res- 
olution, we also performed {dlNs/dV) mapping of an 
optimally doped single crystal at a constant bias volt- 
age of V « (Ad/e) with tip scanning along {001} and 
k || {100}. Periodic conductance modulations consistent 
with the known atomic layer separations in YBCO were 
resolved. Independent imaging on NbSe2 also confirmed 
atomic spatial resolution. 

For the optimally doped and underdoped YBCO, the 
STS exhibited long-range (~ 100 nm) spatial homogene- 
ity and strong directionality, showing a zero-bias conduc- 
tance peak (ZBCP) for k || {110} @, |, |, §|, H, (Fig. 
0(a)); nearly "U-shape" gap features around the zero bias 
for k || {100}, (Fig. 0(b)); and "V-shape" features for 
k || {001}, (Fig. 0(a)). By taking into account a finite 
transverse momentum distribution for the incident quasi- 
particles relative to the normal of the crystalline plane, 
(i.e., by considering a finite "tunneling cone", parameter- 
ized as j3 with a typical value ranging from 15° to 22.5° 
||), the primary features and the directionality of the 
spectra of all optimally doped and underdoped YBCO 
samples were consistent with the generalized Blonder- 
Tinkham-Klapwijk (BTK) theory by Hu 125 1, K ashiwaya 
and Tanaka Q for d x 2_ y 2 pairing [f^ jsj |fj[" 25j, . 

For the Ca-doped YBCO epitaxial films, macroscopic 
spatial variation in the STS at a length scale of ~ 50 nm 
was observed, which correlated with the dimension of the 
growth islands according to images of atomic force mi- 
croscopy j^U , while the STS within each island were spa- 
tially homogeneous, as exemplified in Fig. |l](c). The spa- 
tial modulation of STS may be the result of macroscopic 
phase segregation, similar to the finding in overdoped 
LSCO |l|. Furthermore, the STS of Ca-YBCO exhib- 
ited long-range and symmetric "subgap peaks" , (see Figs. 
0(c) and 0(a)), which differed from the spectral contri- 
butions of local disorder (such as oxygen vacancies), the 
latter generally appeared as short-range and asymmet- 
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FIG. 1: Spatially resolved {dlNs/dV) vs. (V) spectra at 
4.2 K: (a) Optimally doped YBCO. (b) Underdoped YBCO. 
Panels: normalized {100} spectra of an underdoped YBCO 
(upper panel) and (Zn,Mg)-YBCO (lower panel) together 
with a BTK fitting. curve (dashed line), (c) Ca-YBCO. The 
tunneling cone [^, 0| of all spectra ranges from (n/ 12) to 
(•7r/8). For the BTK analysis, each set of the data was normal- 
ized relative to the polynomial fit to the high-voltage back- 
ground conductance, as shown in (b) and (c) by the dashed- 
dotted curve. 



ric "humps" in the spectra. These spectra were in good 
agreement with the (d x i_ y i + s) pairing symmetry ac- 
cording to the generalized BTK analysis |0, ||, 0, |6) , 
with a pairing potential A^ = A^ cos 28k + A s , where 8k 
is the angle of k relative to {100}. We obtained two typ- 
ical sets of spectra with A^ — 17 meV and A s = 9 meV 
for data in Fig. 0(c), and A^ = 13 meV and A s = 6 
meV for the other set of spectra given elsewhere [0. 
On the other hand, d x 2_ y 2 pairing symmetry was con- 
firmed on the optimally doped and underdoped YBCO 
c-axis epitaxial films prepared by two research groups, 
as exemplified in the inset of Fig. 0(a) for a YBCO film 
with T c — 85.0 ±1.0 K. For clarity, we depict the calcu- 
lated c-axis tunneling spectra in Fig. 0(b). Thus, under 
the premise that the spectra are spatially homogeneous 
well beyond the coherence length and mean free path, 
we suggest that the long-range and symmetric subgap 
peaks in overdoped YBCO represent supporting evidence 
for doping-induced variations in the pairing symmetry. 
The substantial s-pairing component (> 20%) may also 
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FIG. 2: (a) Normalized c-axis tunneling spectra of one under- 
doped, one (Zn,Mg)-doped, and one Ca-YBCO, together with 
a BTK fitting curve to the Ca-YBCO spectrum (the thick 
solid line). Inset: c-axis spectrum of an underdoped YBCO 
film, (b) Calculated spectra for different pairing symmetries. 



contribute to the significant enhancement in the critical 
current of bi-crystal Josephson junctions and polycrys- 
talline samples based on the Ca-YBCO |2jq] . 

The maximum value of the d-wave gap, A^, appeared 
to be non-monotonic with the doping level p (Fig. ||(a)), 
whereas the ratio (2Ad/ksT c ) increased with decreasing 
doping, from ~ 7.8 for p w 0.09 to ~ 4.5 for p ss 0.22, as 
illustrated in Fig. ||(b). On the other hand, experiments 
on Bi-2212 found an averaged gap A* that scaled with 
the pseudogap temperature T*, 2 A* ~ 6M B T* @, ||, 
[1G| t \u\, 18, 19|, and that A* increased with decreasing p 
|j,]l|, H|, ||, HI] . We speculate that the superlattice 
and charge-density- wave modulations of the BiO layers in 
Bi-2212 may have complicated the quasiparticle spectra 

@, m, m & Furtherm ° re ' the tendenc y 

of random oxygen distribution in Bi-2212, as opposed 
to the generally ordered oxygen distribution in YBCO, 
may be responsible for the nano-scale variations in A* 
of Bi-2212, ranging from ~ 15 to ~ 70 meV over a few 
nanometers p"T| . 

Although the generalized BTK analysis was suitable 
for deriving the pairing potential and the primary spec- 
tral characteristics, it could not account for the "satellite 
features" associated with many-body interactions [M in 
all YBCO spectra. In Bi-2212, the spectral "dip/hump" 
features were widely observed HJ, ]lf , (l|, |§ [l§] and 
attributed to quasiparticle damping via interactions with 
collective spin fluctuations J27], Eq]. In the strong cou- 
pling limit, the spectral dip is expected to appear at the 
energy u>di P = A + f2 res , where A is the measured gap, 
and £l res is associated with the resonance of collective 
spin excitations J27], |2^| . Defining fl res as the energy dif- 
ference between the primary peak (A) and the dip (cudip) 
in the spectra pa], as indicated by the solid arrows in 
Figs. |l|(b) and]l|(c), we found that £l res in YBCO de- 
creased with decreasing p, as shown in the inset of Fig. 
1(a)- 

In the (Zn,Mg)-YBCO single crystal, microscopic spa- 
tial variations were observed. Two types of impurity scat- 



> 

CD 

< 



< 

CM 



< 



(a) 



> 60 

CD 

E 40 
20 



YBCO 
■(NIS 



Bi-2212 
(SIS) 



ft Krasnov et al 



J 0.0 0.1 0.2 

] Miyakawa et al T ^'-> 
d Renner et al 16 ^~~T"-<&b 



' -FX; 



YBCO 



(Zn.Mg)-YBCO ~ 



Krasnov et al 
M iyakawa et al 1J 
Renner et al 16 




0.00 0.05 0.1 0.1 5 0.20 0.25 0.30 
P 

FIG. 3: (a) Comparison of A d (p) in YBCO with A*(p) in 
Bi-2212. The doping level p is determined from the formula 
T c = T c , max [l - 82.6(p - 0.16) 2 ] § with T c , max = 93.0 K. 
Inset: Comparison of n res (p) for YBCO and Bi-2212. (b) 
Doping-dependent (2A d /k B T c ) for YBCO and (2A*/k B T c ) 
for Bi-2212. The error bar associated with each doping level 
of YBCO covers the range of Ad (or Q r es) obtained from all 
spectra and the uncertainties of the BTK fitting. 



tering spectra were found. One was associated with a res- 
onant scattering at (— 10±2) meV, as shown in Fig. ||, the 
other at (4 ± 2) meV. Assuming the position where a lo- 
cal maximum intensity of a resonant peak occurred as an 
impurity site, we found that the peak persisted over sev- 
eral lattice constants for displacement along either {100} 
or {010}, and the peak intensity decreased rapidly within 
the Fermi wavelength, as shown in the inset of Fig. ^. 
Eventually the usual spectrum was recovered at about 
two coherence lengths (^3 nm) from an impurity site, 
similar to Zn-doped Bi-2212 |^fj| . For displacement along 
other directions, the spectral features became more com- 
plicated j30). Our data have confirmed that the spinless 
impurities such as Zn 2+ and Mg 2+ are strong pair break- 
ers H, H ||, II H, H, i| @- Furthermore, 
the magnitude of the global A^ and Q res in (Zn,Mg)- 
YBCO both became smaller than those of the optimally 
doped YBCO, as shown in Fig. ||(a). The strong suppres- 
sion of fires suggests that the spinless impurity-induced 
moments may be responsible for the suppression of col- 
lective spin excitations [[38] . The single resonance peak at 
the spinless impurity site is also in contrast to the pre- 
dicted double peaks (at ±E ) at the impurity site had 
the T-symmetry been broken due to a complex pairing 
symmetry of either (d x 2_ y 2 + id xy ) or (d x 2_ y 2 + is) p5[ . 
In the context of quantum criticality, our data of doping- 
dependent pairing symmetry is suggestive of a QCP with 
broken C4„-symmetry. However, a small s-pairing com- 
ponent beyond the resolution of our STS might exist in 
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FIG. 4: {001} tunneling spectrum of (Zn,Mg)-YBCO at a 
spinless impurity site and a usual spectrum away from any 
impurity. Two types of impurity scattering were observed: 
one at (— 10±2) meV, the other at (4±2) meV. Inset: Rapidly 
decreasing intensity of the scattering peaks are shown for dis- 
placement along {010} from an impurity site. 



the under- and optimally doped YBCO due to the crys- 
talline orthorhombicity, which would imply that no obvi- 
ous broken symmetry had taken place. For comparison, 
no doping dependence in the pairing symmetry has ever 
been observed in tetragonal Bi-2212. On the other hand, 
the absence of either (d x 2_ y 2 + id xy ) or (d x 2_ y 2 + is) 
pairing does not rule out the possibility of broken T- 
symmetry at a QCP, because certain broken T-symmetry 
states, such as the staggered flux || |3^] and the circu- 
lating current phase pj, cannot be detected directly by 
tunneling spectra. However, to date no universal broken 
symmetry associated with a QCP can be identified for 
all cuprate superconductors. 

In summary, we have demonstrated long-range spa- 
tial homogeneity in the quasiparticle tunneling spectra of 
YBCO samples with a range of doping levels. In contrast, 
STS of the (Zn,Mg)-YBCO exhibit strong pair-breaking 
effects near the Zn and Mg impurities. The quasipar- 
ticle spectral characteristics and the impurity scattering 
effects in optimally doped and underdoped YBCO are 
consistent with predominantly d x 2_ y 2 (> 95%) pairing 
symmetry, whereas those of the overdoped Ca-YBCO ex- 
hibit (d x 2_ y 2 + s) pairing symmetry with a significant 
s-component (> 20%). 
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